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Item  20  (Continued) 

Briefly,  the  cavity  model  theory  developed  previouely  la  refined  so  that 
it  can  be  applied  even  if  Che  dominant  mode  is  not  strongly  excited. 

This  improved  theory  has  been  tested  for  microstrip  antennas  of  two  typical 
geometries:  rectangles  and  circular  disks.  In  both  cases,  the  agreement 
with  the  measured  results  is  found  to  be  remarkable. 

It  is  well-known  that  a  microstrip  antenna  can  produce  circularly  polarized 
(CP)  waves  with  only  a  single  input  post  and  no  phasing  networks.  But,  so 
far,  the  required  dimensions  of  the  antenna  are  determined  by  a  painstaking 
trial-and-error  method,  either  by  experiment  or  with  a  computer.  Since  CP 
operation  is  possible  only  for  a  small  portion  of  the  already  small  bandwidth 
of  the  microstrip  antenna,  the  search  for  those  critical  dimensions  of  the 
aforementioned  method  is  clearly  time-consuming  and  costly.  From  our 
improved  theory,  a  surprisingly  simple  formula  is  derived  to  obtain  the 
necessary  dimensions  of  a  nearly  square  CP  patch  antenna.  The  accuracy  of 
this  formula  is  verified  experimentally  with  almost  perfect  agreement. 

This  work  is  extended  to  elliptical  CP  patch  antennas  and  again  the  same 
close  agreement  with  the  experiment  is  observed. 

ince  thin  microstrip  antennas  are  inherently  narrow  band,  it  is  important 
to  have  a  simple  means  for  tuning  the  antenna  over  a  wide  range.  Based  on 
the  cavity-model  theory,  one  can  expect  that  the  antenna  could  be  tuned 
like  a  cavity.  Some  workers  have  demonstrated  experimentally  the  tuning 
with  a  shorting  post  in  the  patch.  It  is  desirable  to  have  a  theory  which 
can  predict  the  exact  post  location  for  a  given  operating  frequency.  To 
this  end,  our  multiport  theory  is  extended  so  that  a  simple  algorithm  can 
be  applied  to  a  wide  frequency  range,  covering  the  resonant  frequencies  of 
several  modes.  Using  this  extended  theory,  the  operating  frequencies  for 
many  post  locations  are  computed  and  found  to  agree  with  the  measured 
results  within  a  fraction  of  one  percent.  In  fact,  it  is  believed  that  the 
theory  can  be  improved  even  further.  It  is  also  shown  that  CP  can  also  be 
obtained  by  placing  a  post  in  a  square  patch.  Again,  the  theoretical  and 
experimental  results  are  in  almost  perfect  agreement  ..s 

Most  microstrip  antenna  studies  so  far  are  characterised  by  a  simply 
connected  region  in  geometry.  Since  our  theory  implies  that  it  is  mainly 
the  equivalent  magnetic  current  along  the  patch  boundary  that  radiates, 
the  question  has  been  raised  as  to  whether  patches  bounded  by  two  or  more 
closed  curves  would  have  different  radiation  characteristics.  To  obtain 
a  quick  answer  to  this  question,  annulus  and  "quadrulus"  (i.e.,  a 
rectangle  with  a  concentric  rectangular  region  removed)  patches  are 
investigated.  No  significant  difference  is  observed  except,  as  expected,  foi 
the  change  in  resonant  frequencies  and  cross-polarized  component  patterns. 

Previous  investigations  show  chat  the  impedance  locus  of  a  microstrip 
antenna  can  be  varied  over  a  wide  range  simply  by  changing  the  coaxial  feed 
location  in  the  patch,  thus  providing  a  simple  means  for  impedance  matching. 
However,  the  use  of  cable  feed  increases  the  manufacturing  complexity  since 
it  '’an not  be  produced  simplv  by  "printed  circuit"  technique.  To  alleviate 
this  difficulty,  the  method  of  using  a  stripline  feed  through  a  slit  cut 
in  the  patch  is  investigated.  Similar  impedance  variation  is  obtained. 

In  i  microstrip  antenna,  there  is  a  dielectric  loss  in  the  substrate  an  ; 
copper  loss  in  the  cladding.  Manufacturers  usually  provide  die  value  i 
dielectric  loss  tangent  with  no  mention  of  the  copper  loss.  But  in  the 
manufacturing  process,  the  copper-to-dielectric  interface  is  physically 
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Item  20  (Continued) 


and  chemically  treated,  and  as  a  result,  there  is  a  profound  effect  on 
the  copper  loss.  Conventional  methods  of  measurement  cannot  determine 
this  copper  loss  without  altering  the  surface  condition.  Therefore,  a 
new  technique  is  developed  which  can  separate  the  two  losses  in  measure¬ 
ment  without  physically  separating  the  two  materials.  The  result  shows 
that  the  loss  in  copper  cladding  is  substantially  higher  than  what  would 
be  computed  by  using  the  published  conductivity  value  for  pure  copper. 

The  most  challenging  problem  in  microstrip  antennas  is  how  to  broaden 
the  bandwidth.  The  method  of  stacking  two  patches  tuned  to  slightly 
different  frequencies  is  not  very  successful  because  the  small  gain  in 
bandwidth  could  have  been  achieved  much  more  simply  by  using  a  single  patch 
of  double  thickness,  since  both  bandwidth  and  efficiency  increase  with  the 
substrate  thickness.  But  increasing  the  thickness  is  not  desirable  in  many 
applications.  Thus  a  planar  array  of  two  thin  microstrip  elements  of 
slightly  different  dimensions  is  investigated.  Using  the  network  representa¬ 
tion  for  the  elements,  an  efficient  algorithm  for  the  array  is  developed  to 
search  for  an  "optimum"  design  in  the  sense  that  the  bandwidth  is  largest 
for  a  given  SWR  and  simultaneously  maintaining  a  stable  pattern.  Unfortun¬ 
ately,  the  bandwidth  gained  is  not  appreciable.  But  the  algorithm  can  be 
used  to  design  an  array  for  dual- frequency  operation. 

For  a  very  thick  microstrip  antenna,  or  one  operating  at  very  high  frequencies, 
the  cavity  model  may  no  longer  be  valid.  For  this  case,  a  rigorous  theory 
is  formulated  in  terms  of  integral  equations.  Because  of  the  numerical 
complexity,  only  some  preliminary  work  is  performed  up  to  this  point. 
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INTRODUCTION 


This  report  covers  most  of  the  work  performed  under  the  Contract  F19628- 
78-C-0025.  During  the  course  of  the  investigation  activities  in  each  quarter 
have  been  periodically  reported  while  the  accomplishment  of  each  major  task 
has  been  published  in  four  Interim  Reports.  Also  reported  is  a  complete  and 
efficient  computer  program  for  evaluating  all  antenna  characteristics, 
including  input  impedance  locus,  radiation  patterns,  power  lost  and  radiated, 
Q-factor,  and  efficiency,  of  the  widely  used  rectangular  microstrip  antenna 
with  a  single  input  port,  or  two  input  ports,  or  one  input  port  and  the  other 
arbitrarily  loaded.  In  this  final  report,  for  completeness,  a  few  highlights 
(which  may  be  useful  for  the  designers)  in  the  four  Interim  Reports  will  be 
summarized,  while  some  other  findings  which  have  not  been  reported  will  be 
discussed  in  some  detail. 

It  is  well-known  that  an  elliptical  patch  element,  like  a  nearly  square 
patch,  if  designed  and  fed  properly,  can  radiate  circularly  polarized  waves. 
Hitherto,  they  are  determined  only  by  the  painstaking  and  time-consuming  trial 
and  error  method  (12].  Recently,  Shen  [9]  uses  the  elliptical  cylinder  functions 
pertaining  to  the  geometry  for  the  modal  expansion  of  the  field,  and  to  avoid 
the  computational  complication  he  resorts  to  the  approximations  for  those 
functions,  when  the  axial  ratio  is  small.  From  this  a  numerical  search  is  then 
made  for  the  best  axial  ratio  (actually,  in  his  paper  the  best  operating 
frequency  is  sought  for  a  given  axial  ratio).  Clearly,  this  method  of  computer 
searching  is  also  time-consuming  and  costly.  Using  a  perturbation  technique 
in  conjunction  with  the  simple  theory  we  developed  previously,  the  required 
axial  ratio  can  readily  be  determined,  and,  in  fact,  it  is  found  to  be  in  almost 
perfect  agreement  with  the  experiment.  This  is  discussed  in  Section  III.  From 
this  study,  some  general  conclusion  will  be  made. 
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Schaubert,  et.  al  [15]  reported  that  a  microstrip  antenna,  as  a  cavity, 
can  be  tuned  over  a  wide  frequency  band  by  placing  posts  in  the  patch,  and, 
in  fact,  circular  polarization  can  be  obtained  in  a  similar  manner.  But 
their  work  is  mainly  experimental  on  a  trial  and  error  basis.  For  an  analytic 
solution  to  this  problem,  the  multiport  theory  developed  earlier  [ 6  ]  can  be 
applied.  However,  in  all  the  applications  studied  previously,  only  a  narrow 
band  in  the  vicinity  of  the  resonant  frequency  of  any  exciting  mode  is  of  con¬ 
cern,  while  in  this  application,  an  efficient  wide  band  algorithm  is  needed. 

This  is  discussed  in  Section  IV. 

Most  patch  antennas  investigated  so  far  are  bounded  by  a  single 
closed  curve.  Since,  according  to  the  cavity  model  theory,  it  is  the  equivalent 
magnetic  current  along  the  patch  perimeter  that  radiates,  question  has  been 
raised  as  to  whether  patches  bounded  by  two  or  more  closed  curves  would  have 
different  radiation  characteristics.  To  answer  this  question,  patch  antennas  of 
two  different  geometries  are  investigated.  The  first  is  an  annulus  which  can 
be  analyzed  theoretically,  and  the  other  is  a  rectangular  patch  with  a  concentric 
similar  rectangular  region  removed.  The  latter  shall  be,  for  convenience,  called 
"Quadrulus."  The  results  show  that  so  far  as  the  radiation  pattern,  the 
impedance  characteristics,  and  the  bandwidth  are  concerned,  they  are  essentially 
the  same  as  that  of  the  patches  without  holes.  The  most  noticable  difference 
is  the  change  in  resonant  frequency,  and  in  some  cases,  the  magnitude  and  the 
pattern  of  cross-polarized  field  components.  This  is  reported  in  Section  V. 

One  interesting  property  of  the  patch  antenna  is  its  built-in  mechanism 
for  impedance  matching.  As  shown  in  our  previous  investigation,  this  is 
accomplished  by  feeding  the  patch  internally  with  a  coaxial  cable.  But  in  so 
doing,  the  antenna  is  not  strictly  monolithic,  so  that  the  advantage  in  manu¬ 
facture  by  using  the  printed  circuit  technique  is  lost.  For  this  reason  an 

experimental  program  is  launched  to  investigate  whether  the  cable-feed  inside 
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Che  patch  could  be  replaced  by  a  stripline-feed  through  a  slot  cut  in  the 
patch.  Indeed,  similar  impedance  variation  is  obtained.  This  is  discussed 
in  Section  VI. 

The  rigorous  theory  of  microscrip  antenna  has  been  considered  off  and 
on  during  the  last  part  of  the  contract  period.  The  purpose  is  to  prepare  for 
a  major  investigation  in  Che  next  contract  in  which  the  applications  of  micro- 
strip  antennas  at  much  higher  frequency  will  be  studied.  Since  this  problem 
is  closely  parallel  to  that  of  an  infinite  periodic  array  of  patch  elements, 
both  problems  are  examined  simultaneously.  Because  of  the  numerical  complica¬ 
tion,  only  formulation  is  reported  here. 


II.  AN  IMPROVED  THEORY 


Analytic  works  on  microstrip  antennas  can  be  grouped  into  four  different 
approaches.  The  first  due  to  Munson  [1]  and  his  followers  [2]  is  based  on  a 
transmission  line  model,  in  which  two  opposite  sides  of  a  rectangular  patch 
are  considered  as  the  radiating  edges  and  they  are  connected  by  a  low  impedance 
stripline,  namely  the  patch.  They  use  the  impedance  of  a  semi-infinite  parallel- 
plate  guide  derived  approximately  by  Harrington  [3]  as  the  impedance  of  each 
edge.  Clearly  many  questions  could  be  raised  about  this  approach  since  each 
edge  is  not  geometrically  identical  to  the  parallel-plate  guide  in  many  respects. 
Furthermore,  when  a  feed  is  connected  to  one  edge,  as  is  usually  the  case, 
exactly  how  the  impedance  looking  into  the  feed  line  is  related  to  the  edge 
impedance  is  open  to  question.  Obviously,  this  approach  is  difficult,  if  not 
impossible,  to  apply  for  any  general  feed  location  in  the  patch,  and  also 
cannot  be  applied  to  geometries  other  than  rectangle.  It  is  therefore  not 
surprising  to  see  much  disagreement  between  its  predicted  and  measured  impedance 
loci  except  for  some  special  cases. 

A  second  approach  is  purely  numerical  by  modeling  the  antenna  as  a  grid 
of  wires  [4].  While  the  method  seems  to  give  reasonably  good  results,  it  is 
costly  and  provides  little  understanding  of  the  operation  mechanism  of  the 
antenna. 

Taking  advantage  of  the  thinness  of  the  antenna,  a  third  approach  is  made 
by  modeling  the  antenna  as  a  cavity  [5].  It  is  based  on  the  assumption  that 
the  field  structure  under  the  patch  is  not  much  different  from  that  of  a 
corresponding  cavity,  which  is  made  of  the  patch,  the  ground  plane,  and  a 
magnetic  wall  along  the  perimeter.  Once  this  is  accepted,  the  usual  analytic 
procedure  can  be  followed  for  a  solution.  First,  the  field  is  expanded  in  terms  of 
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modes  and  second  from  the  expansion  one  can  find  the  stored  energies,  the 
dielectric  and  copper  losses, as  well  as  the  radiated  power.  From  all  these 
quantities  the  input  impedance  can  be  computed  from  any  one  of  the  three 
formulas:  (a)  admittance  in  terms  of  total  real  power  loss,  stored  energy  and 

input  voltage,  (b)  impedance  in  terms  of  total  real  power  loss,  stored  energy 
and  input  current,  and  (c)  impedance  as  the  ratio  of  input  voltage  to  input 
current.  In  the  theory  developed  previously  [5]  computations  based  on  (a)  have 
given  results  in  close  agreement  with  the  measured.  However  for  the  feed 
location  where  the  dominant  modes  are  not  strongly  excited,  not  only  erroneous 
results  are  obtained,  but  the  computations  based  on  the  above  three  formulas 
give  inconsistent  results.  An  improved  theory  Is  therefore  developed  in  con¬ 
junction  with  a  highly  efficient  numerical  method  [6].  It  yields  accurate 
results  for  any  feed  location  and  also  for  several  critical  applications.  This 
approach  is  simple  and  elegant,  but  also  provides  much  physical  insights  into  the 
mechanism  of  the  antenna  operation,  that  lead  to  some  improved  designs  and 
applications.  However,  this  accomplishment  is  not  made  without  limitations. 

One  of  the  limitations  is  that  it  is  applicable  only  to  a  handful  of  geometries 
which  were  reported  previously  [5].  This  is  not  really  important  since  the 
basic  properties,  such  as  impedance,  patterns,  efficiency,  and  bandwidth,  for 
various  geometries  are  essentially  the  3ame.  The  second,  and  the  most  important, 
limitation  is  that  the  theory  is  useful  only  for  relatively  thin  microstrip 
antennas,  perhaps  on  the  order  of  a  few  hundredths  of  wavelength,  depending 
on  the  accuracy  desired. 

A  fourth  approach  is  the  most  rigorous  one  in  which  integral  equations 
can  be  formulated  by  imposing  the  boundary  conditions  and  the  source  condition. 
Unfortunately,  the  exact  solution  to  these  equations  is  practically  impossible 


to  obtain,  and  even  the  approximate  numerical  solution  is  very  involved.  So 
far  only  some  fundamental  work  has  been  performed  and  hopefully  it  will  lead 
to  a  solution  in  Che  next  phase.  The  advantage  of  this  approach  is  that  it 
will  be  applicable  for  any  antenna  thickness  and  will  also  be  useful  for  the 
design  of  microstrip  surface  wave  antennas.  In  this  Section,  we  shall  summarize 
the  improved  theory  and  leave  the  discussion  of  the  fourth  approach  to  a  later 
section. 


Rectangular  Patch 

First,  consider  the  most  commonly  used  rectangular  microstrip  antenna  as 
shown  in  Fig.  2.1.  The  electric  field  under  the  patch  (assumed  to  be  in  the 
xy-plane  with  dimensions  a  x  b) ,  has  a  z-component  only  and  can  be  represented 
as 
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where  k*  (l-jS^^ )kg,  k^  <*  2irf/v,  f  »  frequency,  v  *  speed  of  light, 

e  =  relative  dielectric  constant  ■  e./e.,  6  cc  =  effective  loss  tangent  =  1/Q 
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cos(mirx/a)  cos(niry/b),  e_  *  1  for  m  *  0  and  2  for  m  i  0,  and  d  is  the  "effective 
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width"  of  a  uniform  strip  of  z  directed  source  current  of  one  amp. 

There  is  an  alternative  representation  based  on  the  modal-matching  to  a 
current  source  of  one  ampere.  It  contains  only  a  single  summation,  and  there¬ 
fore  is  more  suitable  for  numerical  computation.  This  is  shown  below: 
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Circular  Patch 


For  a  disk  patch  of  radius  "a"  and  fed  with  a  uniform  current  source  of 
one  ampere  at  p  -  o’  and  |  +  |  <  y2>  the  electric  fleld  ig  given  below; 
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Radiated  Power  Computation 

Applying  Huygen's  principle  to  the  outer  surface  of  the  cavity  and 
neglecting  the  electric  current  flowing  on  the  outer  surface  of  the  microstrip 
antenna  one  obtains  on  the  magnetic  wall  C  the  Huygen's  magnetic  current  source 

ic  -  2n  x  zE 

z 

where  n  is  the  outward  normal  to  the  H-wall,  and  the  factor  of  2  accounts  for 
the  presence  of  the  ground  plane.  The  electric  vector  potential  of  K  is 


F(r) 


^  -  e0  } 


K(r') 
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e  U  dl(r) 
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and  the  far  field  at  r  is 
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The  total  radiated  power  is 
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Computation  of  Power  Losses  and  Stored  Energies 

The  time-averaged  stored  electric  energy  is  given  by 


W 


e 


(2.7) 


where  t  is  the  thickness,  the  permittivity  of  the  substrate,  S  the  patch 
surface  area,  and  E  is  given  by  Eqs.  (2.2)  and  (2.3)  for  rectangular  and 
circular  patch  respectively.  At  resonance  the  time-averaged  stored  magnetic 
energy  equals  Wg. 
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The  loss  in  the  dielectric  with  loss  tangent  6  is 


me  ,6t 

U 


dS  -  2 u)6  W 


(2.8) 


The  loss  in  the  copper-cladding  with  conductivity  o  is  approximately 


p  rn  — Z —  y  <s  -_Z —  y 

cu  aAty  h  oAty  e 

1/2 

near  resonance,  where  A  »  [2/wya]  is  the  skin  depth,  and  the  factor 
accounts  for  the  surface  area  at  z  ■  0  and  t. 


(2.9) 
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Input  Impedance  Computation 

Define  the  driving  voltage  at  the  feed  as 

V  ■  tE  (at  the  feed,  averaged  over  the  feed  width)  .  (2.10) 

Then  the  antenna  input  impedance  can  be  computed  from  any  one  of  the  following 
three  formulas: 


z  =* 

[V/(-I) ]  at  the  feed 

(2.11) 

z  - 

[P  +  j2u)(Wh  -  We)  ]/ j  ij  2 

(2.12) 

z  - 

|V12/[P  +  j2u»(We  -  Wh)]  , 

(2.13) 

where  P  »  P  +  P,  +  P  .  It  shall  be  noted  that  Eqs.  (2.1)  -  (2.3)  are 
r  a  cu 

derived  for  I  ■  1  amp.  Obviously,  it  is  simpler  to  use  the  first  formula  (2.11) 
for  computing  Z.  Details  for  various  cases  can  be  found  in  the  Interim  Reports 
[6,7]. 
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(b)  Measured  and  computed  impedance  loci. 


6  70  cm 


(INCREASING  FREQUENCY  IS  CLOCKWISE) 


Verification  of  the  Theory 

Two  typical  resultSjOne  for  rectangular  and  the  other  for  circular  micro¬ 
strip  antenna,  are  shown  in  Figs.  2.2  and  2.3.  It  is  seen  that  the  theoretically 
predicted  impedance  characteristics  agree  with  those  measured  remarkably  well  at 
every  frequency  considered. 

Efficiency  Analysis 

Figs.  2.4  and  2.5  show  the  dependence  of  various  power  losses,  including 
the  radiated  power,  on  the  width  "a"  of  a  rectangular  patch  antenna  and  on 
the  substrate  thickness  "t,"  respectively,  when  the  (0,1)  mode  is  excited. 

It  is  interesting  to  see  firstly  that  the  antennp  efficiency  increases  with 
"a"  but  more  slowly  after  a  certain  value, and  secondly,  that  the  efficiency 
increases  sharply  with  t  at  the  beginning,  but  may  even  decrease  for  large  t 
due  to  a  significant  portion  of  the  power  being  converted  into  surface  wave. 

Thus  very  thin  microstrip  antenna  is  not  efficient. 
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Frequency:  I  GHz 
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III .  Circularly  Polarized  Microstrip  Antennas 

.  It  is  well-known  that  circular  polarization  can  be  produced  by  a  microstrip 
antenna  with  a  single  input  port  and  no  phasing  network.  But  all  works  reported 
are  based  on  a  trial-and-error  approach  either  experimentally  [8]  or  with  a 
computer  [9].  Both  of  these  approaches  are  not  only  time-consuming  and  costly, 
but  also  provide  little  physical  insights  into  the  mechanism  of  operation.  In 
this  study  the  conditions  for  generating  circular  polarization  are  formulated 
and  some  surprisingly  simple  but  very  accurate  design  formulas  are  derived. 
Briefly,  the  condtions  are: 

(a)  The  microstrip  antenna  must  be  capable  of  supporting  two  degenerated 
modes;  i.e.,  two  modes  having  the  same  resonant  frequency; 

(b)  the  two  modes  must  radiate  two  orthogonal  field  components  in  the 
direction  normal  to  the  patch; 

(c)  a  proper  perturbation  of  the  patch  geometry  can  be  made  such  that 
the  degenerated  wave  number  will  split  into  two,  say  k^  and  k^,  with 

their  difference  |  k.^  —  i  t0  the  operating  wave  number  k  ratio  equal  to 
the  inverse  of  the  Q-factor  of  the  patch;  i.e., 

|kL  -  k2j/k  =*  1/Q  (3.1) 

(d)  the  feed  must  be  placed  at  an  appropriate  location  for  the  desired 
polarization. 

For  simple  geometries,  such  as  squares  and  circles,  the  above  equation  can 
be  translated  into  the  required  geometrical  perturbation.  There  are  different 


ways  to  accomplish  this,  but  in  the  following  the  energy  perturbation  method 
will  be  presented  because  of  its  simplicity. 
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Perturbation  Theory 

Assume  that  a  magnetic  cavity  of  volume  v  and  resonant  wave  number  k.Q 
is  perturbed  by  removing  (or  adding)  a  small  volume  AV.  Then  the  change  in 
resonant  wave  number  from  k^  to  k^  is  given  by 


kl  ‘  kQ  AW*  “  AWh  -fAv(e!E!2  '  u!Hl2)  dV 

—  - -  *  — - . — ■—  — - -  ■  - 

h  /  ( e 1 E J 2  +  p|H|2)  dv 


W  +  W 
e 


(3.2) 


where  Wg  and  are  time-averaged  stored  electric  and  magnetic  energy  in  v  , 
respectively,  AWg  and  AW^  are  those  in  perturbed  volume  Av,  and  e  and  u  are, 
respectively,  the  permittivity  and  permeability  of  the  medium  in  the  cavity. 

Nearly  Square  Patch 

For  a  rectangular  patch  defined  by  jxj  <  a/z  and  |yj  £  b/2  with 
0  <  a  -  b  *  c  «  a, 
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Since  at  resonance,  W,  =*  W  ,  one  obtains 
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This  implies  the  required  dimension  of  nearly  square  paten  for  CP  is  as 
follows: 

a/b  -  1  +  (1/Q)  .  (3.4) 

This  result  agrees  exactly  with  that  obtained  in  [6]  through  a  different  approach. 
The  surprisingly  high  accuracy  of  this  formula  has  been  verified  experimentally 
as  reported  in  [6]. 

Nearly  Circular  Disk  Patch 

Let  the  semi -major  and  semi -minor  axes  of  an  ellipse  be  a  and  b  respectively, 
and  a  -  b  *  c  <<  a.  Then 

(x/a)2  +  (y/b)2  »  1 

In  terms  of  polar  coordinates  (p,#,z),  the  ellipse  can  be  expressed  approximately 
as 

p  -  [1  -  (c/b)  sin2#]  a  +  0(c2/a2)  .  (3.3) 

Now  consider  first  a  circular  patch  of  radius  a,  excited  for  the  dominant  mode 
(1,1).  The  fields  are 


E2  ”  J1(k^1o)  cos  <j> 

=  -(jk/oiu)  J^k^p)  cos  $ 

Hq  «  -( jk/ajy)  J^tk^p)  sin  <j> 

where  the  modal  wave  number  k^  satisfies  the  characteristic  equation 

J'(kua)  -  0  ,  and  kua  *  1.841  .  (3.6) 

The  change  in  stored  electric  energy  from  a  circular  patch  to  an  elliptical 


patch  is  approximately 
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Similarly,  the  change  in  stored  magnetic  energy  is 
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Thus  when  the  antenna  is  fed  at  x  «  a  (or  -a)  and  y  *  0,  the  (1,1)  mode  will 

be  excited,  and  the  difference  between  the  resonant  wave  number  k  ,  of  the 

el 

elliptical  patch  and  k  of  the  circular  patch  of  radius  a  is ‘given  by 

3 


1  -  3(k  a) 

(c/4a)  - — — j-  «  0.0408  (c/a) 

1  "  (klla) 


(3.8) 


When  the  antenna  is  fed  at  x  *  0  and  y  »  b  (or  -b) ,  the  modal  field 

Ez  *  J^(k^p)  sin<j>  will  be  excited.  For  this  mode,  one  can  either  repeat  the 

above  computation,  or  simply  interchange  a  with  b  and  the  integration  limits 

2 

in  Eq.  (3.7)  to  "from  b  to  b+  (bc/a)  sin  The  difference  between  resonant 

wave  number  k^,  of  the  elliptical  patch  and  of  a  circular  path  of  radius 

b  is  then  given  by 
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ke2  ~  S 


•  -0.0408(c/a) 


From  Eqs.  (3.1),  (3.8),  and  (3.9),  one  obtains 


1  ke2  '  kel  _  «  .  S  "  ka 


k  -  -0.0816(c/a)  +  k 


-0.0816 (c/a)  + 


a  -  b 


0.9184 (c/a) 


(3.9) 


(3.10) 


Hence,  the  required  dimensions  or  axial  ratio  of  the  elliptical  patch  is  given 
by 

c/j  -  1.0887/Q 
or 

a/b  -  1  +  (1.0887/Q)  (3.11) 


Similar  to  a  nearly  square  microstrip  antenna,  when  a  feed  is  placed  along 
the  line  $  -  90°  (or  -90°) ,  LHCP  (or  RHCP)  wave  will  be  radiated  along  the 
2-axis. 

Long  [9]  has  made  measurements  at  around  1.3  GHz  for  two  CP  elliptical 
patch  antennas  with  thickness  ■  0.16  cm  and  t2  **  0.32  cm,  respectively,  and 
both  having  a  ■  4.00  cm.  The  "best"  ratio  (b/a)  for  CP  operation  among  several 
trials  by  him  is  listed  in  Table  1.  Also  listed  are  the  values  predicted  by 
Shen  [9]  who  used  a  computer  to  search  for  the  "best”  frequency  of  operation 
for  a  given  small  axial  ratio. 
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Table  1.  Comparison  of  Measured  and  Predicted  Values  of  b/a  of 
Elliptical  Patch  for  CP 


Measured  Q 
due  to  Long  [10] 


(b/a) 


Exp 

(by  Long) 


Computed 
(by  Shen) 


Computed 
[by  Eq.  (3.11) 


t^  “  0.16  cm 


96.35 


t.,  =*  0.32  cm 


.6.35 


0.984 


0.989 


0.989 


0.976 


0.979 


0.976 


The  values  in  the  last  column  are  computed  using  Long's  measured  Q  which  can 
also  be  computed  readily  from  the  losses  and  the  radiated  power  of  the  circular 
patch  as  discussed  in  Section  2,  or  Ref.  [7].  The  radiation  pattern  and 
impedance  characteristics  of  the  CP  elliptical  patch  antenna  are  very  similar 
to  those  of  the  circular  patch  except  for  a  sharp  cusp  in  the  impedance  locus. 
Therefore,  there  is  no  need  to  use  the  complicated  Mathieu  functions,  nor  their 
approximate  expressions.  It  should  be  noted  that  the  perturbation  method, 
although  yields  results  almost  in  perfect  agreement  with  the  measured  for  both 
cases  discussed  above,  yet  failed  to  predict  the  required  dimensions  accurately 
for  a  square  patch  with  two  opposite  corners  removed  [11].  Apparently,  in 
that  case  new  edge  conditions  and  therefore  new  fringing  field  are  created  as 
the  corners  are  removed.  Therefore,  this  investigation  suggests  that  the 
perturbation  method  is  applicable  whenever  the  edge  condition  is  not  altered 
drastically 

Polarization  Tuning 

In  the  antennas  discussed  above,  a  certain  polarization  can  be  obtained  only 
with  the  feed  placed  at  an  appropriate  location  in  the  patch,  and  therefore  in 
order  to  change  the  polarization  the  feed  must  be  moved  accordingly.  To 
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alleviate  this  inconvenience  the  variation  of  polarizatii.  a  can  also  be  achieved 
for  a  fixed  feed  by  varying  the  DC  biases  of  two  capacicance  diodes  loaded 
at  adjacent  sides  of  a  nearly  square  patch  [11].  If  the  diodes  are  connected 
at  the  nulls  of  the  (0,1)  and  (1,0)  modes,  then  their  effect  on  the  tuning 
will  be  approximately  independent  of  each  other.  An  experimental  sample 
antenna  has  been  fabricated  as  shown  in  the  sketch  in  Figs.  3.1,  in  which  two 
low  pass  filter  circuits  for  each  diode  are  also  shown.  Fig.  3.2  shows  the 
measured  patterns  for  the  LHCP  and  the  RHCP  fields  at  840  GHz. 

The  tuning  can  also  be  accomplished  by  placing  switching  diodes  properly 
located  in  the  patch.  Obviously,  this  technique  is  practical  only  for  very 
few  tuning  steps.  This  will  be  discussed  in  Section  IV. 


Figure  3.2a.  Left  and  right  hand  CP  patterns  of  the  double-tuned  microstrip 
antenna  biased  for  LHCP. 


IV.  WIDE-BAND  MULTIPORT  THEORY 


Schaubert,  Farrar,  Sindoris,  and  Hayes  [15}  have  shown  that  by  loading 
other  ports  of  a  multiport  microstrip  antenna,  the  apparent  resonant  frequency 
can  be  shifted  considerably  away  from  that  of  any  of  the  resonant  frequencies 
of  its  unloaded  counterpart.  Thus,  in  any  modal  expansion  of  the  field,  no 
particular  resonant  mode  will  dominate  at  some  operating  frequencies  in  con¬ 
trast  to  the  case  of  the  unloaded  antenna.  Indeed,  Schaubert,  et.  al.  have 
gxven  experimental  data  indicating  that  by  properly  loading  a  microstrip  antenna 
the  operating  frequency  of  the  radiator  can  be  changed  anywhere  in  a  range  con¬ 
taining  three  different  resonant  frequencies  of  the  corresponding  unloaded  antenna. 
Thus,  a  multiport  theory  which  is  suited  to  such  wide  ranges  of  frequencies  is 
needed  to  adquately  predict  the  behavior  of  loaded  antennas.  Such  a  theory  is 
proposed  in  this  summary. 

This  wide-band  theory  is  a  modification  of  approach  [6].  The  "Z-parameters" 
of  a  multiport  lossy  cavity  which  models  the  antenna  can  be  expressed  as  a  sum 
over  the  resonant  modes  of  the  cavity: 

Z  -  -jknnnt  E  E  TPq/(k2  -  k2  )  ,  (4.1) 

P<1  0  0  m-0  n-0  1,111  1,111 


where 


mn  I  a 


2  r  v 

mn  nn  ] 

.  a  b  i 


Tpq  -  $  (x  ,y  )  <t  (x  ,v  )  jn(mnd/2a) 
mn  ran  p  p  mn  q  q  0 


j q ( x)  *  sin(x) / x  , 


<t  *  (e  e  /ab)  “  cos(mrx/a)  cos(nnv/b) 
mn  om  on 


e  ■ 2  for  m  »  0  and  1  otherwise  , 
on 
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r)Q  ■  377n,  kQ  is  the  free  space  wave  number,  a,  b,  and  t  are  the  patch  dimensions 
and  dielectric  thickness,  respectively,  and  (x  ,y  )  are  the  coordinates  of  the 


In  this  summation,  the  parameter  ”d"  is  the  "effective  feed  width"  which  is 
used  to  model  the  feed  (and  probe  to  which  any  load  may  be  connected  in  the  other 
ports) .  The  current  that  flows  in  the  feed  is  assumed  to  be  a  uniform  thin  band 
of  width  d  of  current  flowing  from  the  ground  plane  to  the  patch.  This  parameter 
has  been  empirically  determined  and  it  appears  to  be  on  the  order  of  five  times 
the  probe  diameter.  According  to  the  theory  reported  in  [6],  the  wave  number 
within  the  cavity,  k,  is  set  equal  to 

1  1/2 

k  =■  erkQ(l  -  j  -) 

where  is  the  refractive  index  of  the  dielectric,  and  Q  is  the  quality  factor 
of  the  antenna  which  can  be  computed  accurately  from  the  formula. 


where 


Qmn  *  +  «  +  A/O 


6  *  dielectric  loss  tangent 
A  *  skin  depth  in  cladding, 


it  k„ak_b  ,  , 

„rad  0  0  1  2 

q  a  - j — _  -  .  r 

ran  e  e  k„t  a 

ora  on  0  mn 


•ir/2  ,2ir  _ 

a  *  sine  de  d$  I F ( 

mn  Jo  Jo 


-F^siniji  +  F^cosi}: 


F  *  (F  cos*  +  F  sin®)  cos8 
ip  x  y 


(k^a)  sinecos$ 

Fx  *  8  - 2 - 2 

(k^asinecos^)  -  (miT) 
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(kQb)  sin6sln$ 


(kQbsin8sin$)  -  (nir)' 


for  the  (m,n)  mode  of  interest.  In  the  limited  band  theory,  reported  in  [6], 
the  inner  sum  of  (4.1)  was  summed  "in  closed  form"  leading  a  one-dimen9ional 
series  whose  convergence  was  further  accelerated  by  a  Rummer's  transformation. 
That  procedure  works  well  as  long  as  the  frequency  band  of  interest  is  about  the 
(m,n)ch  resonant  frequency.  In  the  present  theory,  however,  the  wave  number  k 
in  each  term  of  the  series  in  (4.1)  is  replaced  by  its  own  particular 


k' 

mn 


erk0(1 


V 


1/2 


up  to  a  certain  limiting  wave  number.  This  limiting  wave  number,  k^,  is  twice 

the  wave  number,  k^,  ■  where  f^  is  the  maximum  frequency  at  which  the 

antenna  is  to  be  analyzed  and  c  is  the  speed  of  light. 

Terms  of  (4.1)  with  k  <  k2>  are  summed  directly.  The  remaining  terms 

of  the  series  are  summed  as  follows.  For  k  >  k„  »  2k,  >  2k  *  2t  k„, 

mn  —  z  1  r  o 

i  i  k2 

— =—  -  -~[0.989  +  1.331  At  +  e(k/k  )]  ,  |e(x)  <  0.013 
k*  _  y*  yf-  mn  1  1 

mn  mn  mn 


Therefore,  the  contribution  to  the  infinite  series  in  (4.1)  due  to  m  and  n  such 
that  kmn  _>  k2>  can  be  adequately  approximated  by 


where 


0.939(TQ  -  SQ)  +  1.331  (T  -  S^k4 


T  =  l  TP^/k2J,+2 

1  k  <k,  mn  mn 
mn  2 
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CO  CO 


z  e't 

m*0  n»0 


P9 

mn 


/k 


2V-2 


mn 


t  -  0,1  ; 
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where  '  means  excluding  the  point  (m,n)  ■  (0,0)  from  sum.  Then  the  inner  series 
defining  and  are  summed  in  "closed  form"  and  the  resulting  one-dimensional 
sums  are  evaluated  by  some  simple  approximate  formulas.  The  T's  and  S's  above 
are  constants  dependent  only  on  the  geometrical  and  material  parameters  of  the 
antenna.  Thus, 


Z 

pq 


aQ  -  0.989(TQ  -  SQ)  ,  ^  -  1.331  (Tj_  -  . 


This  single  expression  is  valid  around  every  resonant  frequency  between  D.C. 


and  f  . 


With  the  Z  parameters  thus  computed,  the  driving  point  impedance  of 


loaded  multiport  microstrip  antenna  can  be  evaluated  using  the  known  load  impedance 
matrix  by  the  usual  network  analysis.  Schaubert,  et.  al  [15]  have  made  an 
investigation  of  a  rectangular  patch  with  a  *  9.00  cm,  b  =  6.20  cm,  =  2.55, 
and  t  *  0.16  cm.  The  antenna  is  fed  with  a  cable  at  x'  =  0,  y'  =  1.5  cm  and  a 
shorting  stub  is  placed  at  x  »  a/2  and  y  *  (b/2)  -  l  as  shown  in  Fig.  4.1.  They 
experimentally  determined  the  operating  frequency,  defined  as  one  at  which  the 


SWR  is  minimum  for  each  of  the  stub  positions  l  (normalized  with  respect  to  b/2  ) 
as  shown  in  Fig.  4.2.  Also  shown  is  the  SWR  vs.  stub  position.  Using  the  theory 
and  algorithm  developed  above,  the  operating  frequencies  for  various  stub 
positions  can  be  predicted  as  shown  in  Fig.  4.3.  Comparing  this  with  that  in 
Fig.  4.2,  one  will  find  an  excellent  agreement  except  for  (1)  a  small  shift, 
about  0.847TS,  in  the  frequency  axis,  and  (2)  a  small  discrepancy  near  d  =  b/2. 

It  is  believed  that  they  are  caused  by  the  coraplicaL=i  fringing  effect  around 
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OPERATING  FREQUENCY  (MHz) 


1400  I - 1 - 1 - 1 - 1 1 

0.00  0.25  0.50  0.75  1.00 

POST  LOCATION  (2  l/b) 


Figure  4.2.  Plot  of  measured  operating  frequency  vs.  post  position. 


SWR 


Theoretically  predicted  operating  frequency  va.  post  position 


the  feed  and  the  shorting  stub  when  they  are  near  the  edge  of  the  patch.  They 
could  be  accounted  for  by  using  appropriate  values  for  the  effective  width. 

Fig.  4.4  shows  the  measured  and  computed  impedance  loci  for  a  square 
patch  for  CP  operation.  The  agreement  is  nearly  perfect.  It  may  be  noted 
that  in  this  case  both  the  feed  and  shorting  stub  are  inside  the  patch. 


V.  ANNULUS  AND  QUADRULUS  MICROSTRIP  ANTENNAS 

Most,  if  not  all,  microstrip  antennas  investigated  extensively  are  char¬ 
acterized  in  geometry  by  that  the  patch  is  a  simply  connected  region.  In 
other  words,  the  patch  is  bounded  by  a  single  closed  curve.  Since  it  is 
generally  understood  that  it  is  the  equivalent  magnetic  current  along  the 
boundary  curve  that  radiates,  one  may  speculate  whether  a  patch  bounded  by 
two  or  more  closed  curves  would  have  fundamentally  different  radiation  patterns. 
To  obtain  a  quick  answer  to  this  question,  patches  of  two  simple  geometries  are 
investigated;  namely,  an  annulus  and  a  "window"  shaped  rectangle  (i.e.,  a 
rectangle  with  a  concentric  rectangular  region  removed.  For  convenience,  it  will 
be  called  "quadrulus") .  The  former  can  be  studied  analytically  [5]  while  for 
the  latter  only  numerical  method  is  possible.  As  will  be  seen  later  from  the 
measured  results,  their  impedance  loci  and  radiation  patterns  are  not  basically 
different  from  those  of  the  patches  of  simply  connected  regions.  Therefore,  it 
seems  Chat  a  major  effort  for  the  theoretical  analysis  of  these  antennas  at 
this  moment  is  not  justified. 

Figs.  5. 1-5. 3  show  the  measured  impedance  locus  and  radiation  patterns 
of  a  circular  disk  patch  antenna  with  radius  equal  to  6.75  cm  and  fed  at  the 
edge  of  the  disk.  They  should  be  compared  with  the  corresponding  characteristics 
shown  in  Figs.  5.4  -5.6  for  an  annulus  patch  of  radii  6.63  cm  and  1.77  cm,  also 
fed  at  the  edge  of  the  outer  circle. 

Fig.  5.7  shows  the  geometry  of  the  rectangular  patch  antenna  and  the 
"quadrulus"  patch  antenna  investigated.  The  impedance  loci  and  patterns  for 
several  lowest  order  modes  of  both  patches  are  shown  ir.  Figs.  5.8-5.17.  It 
should  be  noted  that  the  power  level  was  not  maintained  constant  through  the 
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Figure  5.1  Impedance  locus  of  a  circular  disk  microscrip  antenna 
6.75  cm,  fed  at  the  edge  of  the  disk. 
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Figure  5.2  Radiation  pattern  in  the  $  *  0a  plane  of  the  disk  antenna  stated 

In  Fig.  5.1  at  798  MHz  (Dotted  pattern  is  for  the  cross-polar ized 
component) . 
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Impedance  loci  of  an  annulus  microstrip  antenna  with  radii  equal 
to  6.63  cm  and  1.77  cm  fed  with  a  cable  at  the  outer  edge. 


measurements.  All  antennas  are  made  of  Rexolite  2200,  I/Lo"  thick.  It  is 
seen  that  qualitatively  they  all  behave  nearly  the  same,  and  no  one  has  a 
significantly  wider  bandwidth.  On  the  other  hand,  this  finding  suggests  some 
design  flexibility  for  the  case  where  a  certain  shape,  for  example  rectangle, 
can  not  be  made  to  fit  a  vehicle.  However,  Che  question  of  how  the  resonant 
frequency  is  related  to  the  dimensions  of  the  patch  of  a  given  shape  can  not 
always  be  easily  answered. 
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Impedance  loci  of  the  rectangular  microstrip  antenna  shown  in 
Fig.  5.7  for  the  first  four  modes. 
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Figure  5.9.  Radiation  pattern  In  the  *  -  0*  plane  of  the  rectangular  micro- 
strip  antenna  shown  in  Fig.  5.7,  at  612  MHz  (dotted  pattern  is 
for  the  cross-polarized  component) . 


Figure  5.11.  Radiation  pattern  In  the  p  •  0”  plane  of  the  rectangular 
microatcip  ancenna  shown  in  rig.  5.7  ac  1.110  GHz  (dotted 
pattern  is  for  cross-polartzei  component). 


Figure  5.14.  Radiation  pattern  in  the  $  »  0'  plane  of  che  quadrulus  micro 
strip  antenna  ehovn  in  Fig.  5.7  at  546  MHz  (dotted  pattern  l 
for  cross-polarized  conponent) . 
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Figure  5.16.  Radiation  pattarn  In  the  «  *  0'  plane  of  the  quadrulus  micro- 

strip  antenna  shown  In  Fig,  5.7  at  1.080  GHz  (dotted  pattern  Is 
for  the  croaa-polarlzed  component). 


VI.  STRIPLINE  FEED  TECHNIQUES 

Microstrip  antennas  can  be  fed  with  a  coaxial  cable,  a  stripline,  or 
j  even  with  a  device  of  electric  of  magnetic  induction.  The  advantages  of 

I 

i  the  cable  feeding  are  that  the  feed  can  be  located  anywhere  in  the  patch 

(Fig.  6.1)  to  obtain  a  desired  impedance  characteristic  and  that  the  cable  can 
be  placed  under  the  ground  plane  so  that  the  coupling  between  the  feed  and  the 
antenna  patch  is  eliminated.  But  ths  disadvantage  is  that  the  structure,  in 
contrast  to  the  stripline  feed,  is  not  completely  monolithic  and  becomes  more 
difficult  to  produce.  For  this  reason,  an  experimental  investigation  has  been 
made  to  determine  the  feasability  of  using  a  stripline  to  feed  at  an  interior 
point  of  c:he  patch  through  a  slot  cut  In  the  patch  as  shown  in  Fig.  6.1  (bottom) 
for  three  different  positions.  The  impedance  loci  of  this  antenna  are  shown 
in  Figs.  6.3  -  6.5  for  various  dimensions.  For  comparison  the  impedance  locii 
of  the  patch  fed  with  a  cable  at  the  same  three  different  positions  are  shown 
in  Fig.  6.2.  Although  they  are  not  identical,  the  possibility  of  using  stripline 
feeth  for  obtaining  various  impedance  locus  is  obvious. 


Impedance  loci  of  Che  rectangular  microstrip  antenna  fed  with 
a  cable  at  three  different  locations:  (1)  y’  -  3.11  cm; 

(2)  y'  -  2.29  cm;  and  (3)  y'  -  0.76  cm. 
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Figure  6.3.  Impedance  loci  of  a  rectangular  microstrip  antenna  fed  with 
a  stripline  at  three  different  depths:  (1)  D  »  3.11  cm; 

(2)  D  -  2.29  cm;  and  (3)  D  -  0.76  cm,  with  gap  G  -  0.254  cm 
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Figure  6.4.  Impedance  loci  of  a  rectangular  microstrip  antenna  fed  with 
a  stripline  at  two  different  depths:  (2)  D  -  2.29  cm  and 
(3)  D  -  0.76  cm,  with  Gap  G  »  0.559  cm.  The  locus  for  depth 
(1)  D  •  3.11  cm  is  too  close  to  the  jrj  =*  1  circle  to  be 
useful. 
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VII.  EVALUATION  OF  LOSSES  IN  MICROSTRIP  ANTENNA  MATERIALS 


Double  copper  cladded  laminates  are  the  most  widely  used  material  for 
microstrip  antennas.  In  this  material  there  are  losses  in  the  dielectric  and 
losses  in  the  copper,  usually  specified  by  the  loss  tangent  and  the  conductivity, 
respectively.  Commonly  the  manufacturers  supply  only  the  data  of  loss  tangent 
with  no  mention  of  the  copper  loss.  A  naive  user  may  consult  a  handbook  for 
the  copper  conductivity  for  its  loss,  not  knowing  that  the  treatment  of  the 
copper  and  dielectric  interface  in  the  manufacturing  process  could  change  the 
effective  copper  loss  drastically.  Unfortunately,  conventional  methods  of 
measurement  cannot  determine  this  loss  without  altering  the  surface  condition. 
Therefore,  a  new  technique  is  needed  so  that  the  two  losses  can  be  separated 
without  physically  separating  the  copper  from  the  dielectric  substrate. 

Making  use  of  the  fact  that  the  dielectric  loss  depends  on  the  volume , 
thus  the  thickness  of  the  substrate,  while  the  copper  loss  on  the  cladding 
surface  area,  one  can  obtain  the  following  simple  relation: 


where  Q  is  the  qualifying  factor  of  a  bona-fide  cavity,  made  of  a  piece  of 
the  material  with  all  the  edges  closed  with  copper  strips;  t  and  6  are  the 
thickness  and  loss  tangent  of  the  substrate,  respective;  and  A  is  the  effective 
skin-depth  of  the  copper  cladding.  The  above  equation  describes  a  straight  line 
in  the  A5-plane  with  slope  -1/t  and  intercept  1/Q  on  the  6-axis.  By  constructing 
two  or  more  cavities  from  the  same  material  (preferrably  from  the  same  batch) 
but  of  different  thickness  and  then  measuring  their  Q-factors  and  thicknesses, 
one  can  solve  for  the  two  loss  parameters  6  and  6.  Various  values  of  thickness 
can  be  obtained  by  stacking  and  clamping  together  two  or  more  layers  of  the 
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material  cut  out  from  the  same  sheet  with  the  interlayer  copper  claddings 
removed.  In  our  measurement,  it  is  observed  that  the  value  of  Q-factor  may 
vary  over  a  small  range  when  measurements  are  repeated  day  and  night  over 
a  period  of  time,  but  does  not  vary  significantly  when  the  small  air  gap 
between  layers  are  changed.  The  variation  could  be  caused  by  the  change  in 
physical  conditions,  such  at.  temperature,  moisture,  most  probably  in  our  case 
the  equipment  instability,  etc.  Because  of  these  causes  and  also  because  of 
possible  inhomogenieties  in  the  material,  many  measurements  may  not  yield  the 
same  values  of  A  and  5.  Thus,  it  can  only  be  assumed  that  the  centroid  of 
the  commonly  intersected  region  as  shown  in  Fig.  7.1  would  be  likely  to  give 
the  true  values.  Table  2  shows  the  results  of  an  extensive  measurement  for 
Rexolite  2200  at  various  frequencies.  The  conductivity  calculated  from  them 
is  about  3.64  x  107  mho/M,  which  is  only  5/8  of  the  published  conductivity  of 
pure  copper,  namely  5.8  x  107  mho/M.  It  is  found  that  by  using  this  value  the 
input  impedances  of  many  microstrip  antennas, made  of  Rexolite  2200  and  computed 
from  the  theory  in  Section  II,  agree  with  the  measured  almost  exactly.  This 
perhaps  can  be  regarded  as  an  independent  verification  of  the  method. 
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Figure  7.1.  Loss  tangent  vs.  skin  depch  for  Rexolite  2200  with  substrate 
thicknesses  of  1/16,  1/8,  and  3/16  inch. 


Loss  tangent  6  and  skin  depth  L  of  Rexclite 
2200  at  various  frequencies. 
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.00111 
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.000218 

VIII.  ARRAY  STUDY 


A  major  limitation  of  the  microstrip  antenna  is  narrow  bandwidth,  a 
characteristic  of  all  resonant  structure.  Several  attempts  for  broadening 
the  bandwidth  have  been  made  by  many  workers  [13,14]  with  little  success.  Our 
early  trials  on  log-periodic  arrays  of  microstrip  elements  were  not  successful 
because  first  it  is  impossible  to  achieve  a  criss-cross  feed  for  grounded 
elements  in  order  to  fulfill  the  backfire  conditions,  and  second  an  antenna 
does  not  necessarily  have  frequency-independent  performance  even  if  its  structure 
satisfies  the  so-called  scaling  principle.  The  method  of  stacking  two  elements 
of  slightly  different  dimensions  have  been  reported  [13].  A  close  examination 
of  this  work  shows  that  the  small  increase  in  bandwidth  so  gained  could  have 
been  realized  more  easily  by  using  thicker  substrate  for  a  single  element 
[see  Section  II,  Fig.  2.5]  without  the  complicated  feeding  problem  as  in 
the  stacking  arrangement.  The  obvious  drawback  of  this  approach  is  that  the 
antenna  is  no  longer  thin  and  may  protrude  out  of  the  ground  plane  with 
undesirable  height.  (In  our  previous  report,  it  was  shown  that  the  antenna 
could  be  mounted  flushly  with  the  external  ground  plane.  But  in  this  arrangement, 
the  antenna  would  protrude  inward.)  An  alternative  approach  consists  of  a  small 
array  of  two  elements  of  slightly  different  dimensions.  This  is  inspired  by  the 
multiple-tuning  circuit  for  a  broadband  operation  and  has  been  experimental ly 
explored  by  Pues,  Vandensande  and  van  de  Capelle  [14].  The  objective  of  our 
study  is  to  make  use  of  the  theory  discussed  in  Section  II  to  systematically 
develop  a  design  algorithm  such  that  the  best  design  will  be  obtained  without 
the  painstaking  trial-and-error  method  to  actually  test  many  antennas. 

A  computer  analysis  is  conducted,  modeling  the  array  elements  by  their 
equivalent  networks  and  using  the  transmission  line  theory  to  calculate  the 


input  impedance  at  the  array  feed.  From  this  the  currents  at  each  element 
can  be  determined  and  then  the  radiation  pattern  of  the  array  is  calculated. 
The  design  goals  are  a  stable  pattern,  good  efficiency,  and  low  SWR  over  a 
wider  bandwidth  than  possible  with  a  single  element.  In  the  following  the 
network  representation  for  the  antenna  element  will  be  discussed  first. 


Network  Representation  for  a  Microstrip  Antenna  Element 

The  input  impedance  of  a  microstrip  antenna  element  can  be  represented 
conveniently  in  a  Foster  expansion.  To  illustrate  this,  consider  a  rectangular 
patch  discussed  in  Section  II,  where  it  was  shown  that  for  input  current  of 
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Then  Equation  (8.1)  can  be  written  in  a  more  interesting  form: 
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As  microstrip  antennas  are  narrow  band  devices,  for  operation  in  the 

frequency  band  of,  say,  mode  (M,N),  G  («)  can  be  replaced  by  G  (w„  .,). 

mn  nm  m,  n 

In  doing  so,  Eq.  (8.6)  represents  a  Foster  expansion  of  the  input  impedance 
function  of  a  network  consisting  of  an  infinite  number  of  parallel  resonant 
circuit  in  series.  Furthermore,  since  Lmn  decreases  with  increasing  modal 
indices,  the  infinite  number  of  high  order  Foster  sections  can  simply  be  com¬ 
bined  to  form  a  single  inductance  in  series  with  the  remaining  circuit  elements 
as  shown  in  Fig.  8.1a. 

In  case  that  the  antenna  is  operating  near  a  resonant  frequency  <a_,  which 

MN 

is  well  separated  from  all  other  resonant  frequencies,  a  further  simplification 
of  the  network  representation  is  possible,  as  shown  below: 


Z  a 


WV  +  JuCMN  “  j/wLMN 


+  jwL' 


(8.7) 


where 


L' 


mn 


m,n^M,N  ^  -  a^N 
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This  is  simply  the  impedance  of  a  parallel  resonant  circuit  in  series  with  an 
inductance  L'  as  shown  in  Fig.  8.1b,  and  this  also  explains  why  the  impedance 
locus  is  shifted  to  the  inductive  region. 

Two  Element  Array 

Two  rectangular  patch  elements  of  slightly  different  dimensions  are 
considered.  As  reported  previously  [5] ,  the  mutual  coupling  effect  between 
two  elements  is  not  significant  unless  they  are  unusually  close.  Thus, 
in  general  the  array  design  is  essentially  a  circuit  problem  since  each  element 
can  be  adequately  represented  by  a  circuit  in  Fie.  8.1b  with  impedance  given 
by  Ea.  (8.7).  What  requires  is  mainlv  the  feed  network  design.  To  keep  the 
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Figure  8.1  Circuit  representations  for  nlcrostrip  antenna. 

(a)  General  network  representation  for  frequency  band  about 
a  resonant  mode. 

(b)  A  simplified  representation  for  the  band  about  the  resonant 
frequency  of  (M.N)th  mode  when  it  is  well  separated  from 
all  others. 


structure  monolithic,  simple  stripline  feed  is  considered  as  shown  in  Fig.  8.2 
The  relative  excitation  currents  to  the  elements  as  well  as  the  input  impedance 
at  the  array  feed  point  can  thus  be  expressed  in  terms  of  the  dimensions  of 
each  patch  (a  and  b) ,  feed  location  to  each  element,  and  the  feed  network 
parameters,  namely  the  line  length  d  and  line  width  W  for  each  of  three  line 
sections.  Once  the  currents  are  found  the  radiation  pattern  and  the  efficiency 
of  the  array  can  be  computed.  A  computer  program  is  written  to  evaluate  all 
these  quantities.  The  following  few  parameters  are  chosen  a  priori:  the 
substrate  thickness  t  of  the  material  is  0.15  cm;  the  line  width  is  determined 
for  a  50-ohm  characteristic  impedance;  the  patch  dimensions  are  chosen  to 
resonate  at  two  adjacent  frequencies,  for  example,  795  MHz  and  805  MHz.  By 
observing  the  impedance  locus  and  radiation  pattern  as  the  line  lengths 
d^\  d^,  and  d^  are  varied,  hopefully  we  can  obtain  an  optimum  design. 

Before  searching  for  the  optimum  design,  the  theory  as  well  as  the  computer 
program  are  checked  against  the  experiment.  A  typical  result  is  shown  in 
Fig.  8.3.  The  agreement  is  only  fair.  Considerable  effort  has  been  expended 
to  determine  the  possible  causes  of  the  discrepancy.  This  is  found  partly  due 
to  the  discontinuities  of  the  striplines  and  the  stripline-to-coax  connector 
at  the  feed  point,  and  mostly  due  to  the  inaccuracy  of  our  equipment  which 
is  over  fifteen  years  old.  Inhomogeneity  in  the  material,  and  in  accuracy 
in  the  construction  are  also  possible  contributing  factors.  However,  the 
computer  program  seems  to  be  capable  of  at  least  predicting  the  general 
behavior  of  the  array. 

The  computer  program  is  finally  used  to  search  for  the  best  design. 

Unfortunately,  it  is  found  that  Che  case  for  a  broadband  in  impedance 

characteristics  doe  not  give  a  stable  pattern  over  the  band  (i.e.,  the  beam 

may  tilt  from  one  side  to  the  other).  It  is  a  simple  matter  to  use  our  program 

to  evaluate  the  two  element  array  reported  by  Pues,  Vandersande,  and  van  de 

Capelle  (14].  The  result  shows  that  the  beam  is  skewed  from  the  broadside  bv 
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Figure  8. 
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MEASURED  LOCUS 


cf'*  ■  90,  Q(2)  ■  91,  0.  «  23-6  cm, 

(I)  (2)  (3) 

d  »  49*8  cm,  d  ■  64-5  cm,  d  •  9*31  cm, 

0)  (2)  (3) 

w  »  w  ■  w  ■  0*44  cm,  t  ■  0*15  cm. 


l  «  0*0036  cm 


3.  Measured  and  computed  Impedance  loci  of  the  two-element  array. 
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a  considerable  amount.  Therefore,  cheir  claim  of  wide  bandwidth  is  not  sub¬ 
stantiated.  Although  our  attempt  to  broaden  the  bandwidth  of  a  two-element 
array  is  not  completely  successful,  it  is  believed  the  computer  program  could 
be  used  for  the  design  of  an  array  for  dual,  or  even  multiple  frequency 
operation.  It  is  further  believed  that  the  program  can  be  expanded  to  study 
symmetrical  arrays  of  three  or  more  elements  in  which  case  a  stable  pattern 
could  be  maintained. 


IX.  GENERAL  THEORY  AND  INFINITE  ARRAY 

Early  interest  of  the  microstrip  antennas  has  been  focused  on  their  com¬ 
pactness  in  structure,  and  the  investigations  are  therefore  mostly  confined 
to  thin  substrates.  Taking  advantage  of  this  special  property  one  is  able 
to  develop  a  simple  theory  based  on  the  cavity  model  which  turns  out  to  be 
surprisingly  accurate.  As  in  the  development  of  all  scientific  studies, 
progress  and  inventions  are  often  made  by  extending  one's  knowledge  to  new 
situations.  It  is  therefore  natural  to  expand  our  investigation  to  "thick" 
microstrip  antennas,  or  to  microstrip  antennas  to  be  operated  at  much  higher 
frequencies.  From  our  investigation  this  appears  to  be  the  only  effective  way 
to  broaden  the  bandwidth.  It  seems  that  this  purpose  alone  justifies  a  major 
effort  to  develope  a  general  theory  for  microstrip  antennas. 

If  the  substrate  is  not  sufficiently  thin,  the  assumption  of  vanishing 
transverse  E-components  and  vanishing  longitudinal  H-component  with  respect  to 
the  normal  of  the  patch  is  no  longer  valid,  and  it  follows  that  the  cavity 
model  with  a  magnetic  wall  along  the  patch  perimeter  would  fail.  One  is  there¬ 
fore  forced  to  formulate  the  problem  on  the  satisfaction  of  the  following 
conditions:  (a)  the  boundary  conditions  on  the  patch  and  ground  plane;  (b)  the 

field  continuity  condition  over  the  air-dielecf ric  interface,  and  (c)  the  source 
condition  at  the  feed.  The  last  condition  is  often  overlooked  by  many  workers. 
In  formulating  the  problem  in  this  manner  one  finds  that  it  is  closely 
related  to  that  of  an  infinite  periodic  array  of  microstrip  elements.  As 
stated  before,  mutual  coupling  between  microstrip  elements,  at  least  for 
moderately  thin  substrate,  is  generally  not  strong;  therefore,  the  input 
Impedance  of  any  element  in  the  array  will  approach  that  of  a  single  isolated 
element  as  element  spacing  increases.  This  provides  another  approach  to  the 
single  element  problem;  in  fact,  numerically  speaking,  this  approach  appears 
to  be  somewhat  simpler. 


70 


Geometry  of  a  rectangular  microstrip  antenna  fed  with 
stripline  at y  ■  0,  x  -  c  to  d. 


During  Che  course  of  this  contract,  the  cwo  problems  just  stated  have 
been  studied  off  and  on.  Because  of  the  numerical  complication,  no  quick 
solution  is  anticipated.  The  purpose  is  to  make  some  preliminary  study  which 
hopefully  will  lead  to  a  more  vigorous  investigation  in  the  next  contract  period. 


Simple  Element 

We  shall  consider  a  rectangular  element  of  dimensions  a  x  b.  There  are 

several  ways  to  incorporate  the  source  condition  into  the  formulation.  At 

present  we  assume  a  stripline  feed  at  the  y  ■  0  edge  with  an  input  current  of 

1  amp  distributed  over  the  line  width  x  -  c  to  x  -  d  as  shown  in  Fig.  9.1. 

Let  the  current  density  on  the  patch  be  xf  (x,y)  +  yf  (x,v)  and  its  Fourier 

x  y 

transform  be  xF  (p,q)  +  $F  (p,q).  Then  the  following  set  of  integral  equations 
x  y 

for  Fx(p,q)  and  F^(p,q)  can  be  obtained: 
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This  set  of  integral  equations  can  only  be  solved  numerically.  An  initial 
attempt  has  been  made,  but  the  computation  appears  to  be  so  involved  that  a 
major  effort  is  needed. 

Infinite  Periodic  Array 

Consider  a  periodic  array  of  infinite  number  of  rectangular  patch  elements, 
each  having  the  same  dimensions  and  excitations  as  before,  and  all  arranged  in 
rectangular  lattices  with  periodicity  A  and  3  along  the  x-and  y-axis, 
respectively.  Assume  that  all  elements  are  fed  in  phase  and  the  current  in 
any  one  element  to  be  xJx(x,y)  +  yj^(x,v).  Since  the  current  must  also  be 
periodic  with  period  A  end  3  along  the  x-  and  y-axis,  it  can  be  expanded  in 
Fourier  series  with  coefficients: 


where 


Fx(m,n) 

y 


j{  J  (x, y)  e"3*"**3""7  dxdv 
Ax  3  X 

V 


K  =  2im/ k 
m 


n  =*  2nn/B 
n 

b  <  B 
,'3 


a  <  A 


Since  Jy(x,0)  »  l/(d  -  c)  for  c  <  x  <  d. 


Fy(m) 


d  -j{x 
e  dx 


x-c 


-j£  d  -U  d 
(e  m  -  e  m  )/£ 


After  applying  the  boundary  conditions,  the  following  two  equations  for 

F  (m,n)  and  F  (m,n)  are  obtained: 
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where  all  G-functions  are  the  same  as  before  if  p  and  q  are  replaced  by  ^ 
and  n^,  respectively.  Comparing  the  above  set  of  equations  with  that  for  the 
simple  element,  one  notices  the  close  similarity.  But  the  summations  being 
denumerable,  from  numerical  point  of  view,  it  seems  to  be  simpler  to  consider 
the  infinite  array  problem  first.  It  could  be  solved  approximately  by 
enforcing  the  equality  of  these  equations  at  a  sufficient  number  of  points 
(x.,y.)  for  the  equal  number  of  the  truncated  unknowns  F^fm.n)  and  F^Cm.n). 
Alternatively,  one  can  assume  Jx(x,y)  and  J^(x,y)  be  expressed  in  terms  of  a 
set  of  judiciously  chosen  basis  functions  and  then  solve  for  their  expansion 
coefficients.  Some  simple  basis  functions  have  been  investigated  but  are  found 
nor  suitable  because  of  the  convergence  problem.  But  this  initial  effort  has 
Provided  us  with  much  information  and  experience  on  the  numerical  aspect  of 
•  :-.e  problem.  This  investigation  should  certainly  be  continued. 
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